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ABSTRACT 



Context. Recent studies on stellar evolution have shown that the properties of compact objects strongly depend on metallicity of the 
environment in which they were formed. 

Aims. In this work, we study how the metallicity of the stellar population can affect unresolved gravitational waves background from 
extragalactic compact binaries. 

Methods. We obtain a suit of models using population synthesis code, estimate the gravitational wave background they produce and 
discuss its detectability with second (Advanced LIGO, Advanced Virgo) and third (Einstein Telescope) generation detectors. 
Results. Our results show that the background is dominated by binary black holes for all considered models in the frequency range 
of terrestrial detectors, and that it could be detected in most cases by Advanced LIGO/Virgo, and with Einstein Telescope with a very 
high signal to noise ratio. 
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1. Introduction 

The coalescence of two neutron stars (BNS), two black holes 
(BBH) or a neutron star and a black hole (NSBH), are the 
most promising sources of gravitational waves (GWs) for ter- 
restrial detectors, due to the huge amount of energy emitted 
in the last phase of their inspiral trajectory, the merger and 
the ringdown. The coalescence occurs after two massive stars 
, in a binary system have burned all their nuclear fuel, have 
evolved into red giants, and the cores have collapsed, possi- 
bly after supernova explosions, forming a bound system of two 
compact objects (neutron stars or black holes) inspiralling each 
other due to the emission of GWs. The number of massive sys- 
tems that remain bounded after two supernova explosions (or 
prompt core-collapse) is uncertain, as well as the time to coales- 
cence (the delay), or intrinsic parameters such as masses, spins, 
eccentricity, which depend on complicated evolution scenario 
involving common envelope and mass transfer. The distance 
probed with actual interferometers is about 30 Mpc for BNSs, 
but the next generation of detectors such as Ad LIGO/Virgo 
Q should be taking data with a sensitivity approximately 10 
times greater, pushing the horizon up to about 450 Mpc. With 
the third generation interferometer Einstein Telescope (ET) 
dPunturo et al.l 1201 Ol) . the horizon of compact binaries is ex- 
pected to reach cosmological distances of z > 2 - 3, where it 
may become possible to study the evolution of the sources over 
redshift. Besides the emission produced by the coalescence of 
the nearest binary systems, the superposition of a large num- 
ber of unresolved sources at high r edshifts will prod u ce a back- 
ground of gravitational wave s (see Regimbap (l2011l);IZhu et al.l 
(1201 11) : iMarassi et al.l d20lTl) : iRosadol (boi lh : IWu et al.l (120121) : 



Regimba u et all (120121) : IZhu et all (1201 2l) . for the most recent 
studies) that may dominate over the cosmological background in 
the range 10 - 1000 Hz where terrestrial detectors are the most 
sensitive. In the Laser Interfe rometer Space Antenna (LISA) 
band between 10~ 4 - 0.1 Hz dBenderfll998h . it is expected to 
remain after substraction of the galactic foreground from white 
dwarf binaries. 



Recent studies have shown that the properties of the compact 
object binaries strongly depend on the metallicity of the stellar 
population. The formation rate of binaries containing black holes 
sharply increases with decreasing metallicity, as indicated bot h 
by observations dBulik & Belczvnskil|2010t iBulik et al l 1201 ll) . 
and by binary population synthesis dBelczvnski et a l. 2010a). 
Moreover, it has been shown th at at low metallicity the typical 
mass of a black holes increases dBelczvnski et alfeOlObl) . While 
for the solar metallicity the maximum mass of a black hole is 
about 10-15 M Q , for metallicity of about 10% solar value it rises 
to above 30 M Q . The large redshift population of compact object 
merger originates in populations of stars with various metallic - 
ities including very low ones. Thus investigation of the gravi- 
tational wave background must take into account the effects of 
metallicity. 



1 http://www.ligo.caltech.edu/docs/IVl/M060056-10.pdf, 
https://wwwcascina.virgo.infn.it/advirgo/docs.html 



In this paper, we present a study of the effect of the metallic- 
ity on the GW background from compact binaries (BNS, NSBH 
or BBH). In section II we describe how we use the simulation 
code Star/Track to derive sets of intrinsic parameters, in section 
III we discuss the calculation of the coalescence rate as a func- 
tion of redshift, in section IV we derive the GW background and 
describe the numerical simulations, in section V we present our 
results, and finally in VI we summarize our main conclusions. 
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2. Simulations of a population of compact binaries 
with Star Track 

For a simulated population of compact binaries, the intrinsic 
parameters, such as the masses mi and 1112 of the two compo- 
nents, the eccentricity e, the delay td are selected using the bi- 
nary evolution code StarTrack. The masses of the compact ob- 
jects are determined by the masses of the two progenitors and 
by mass loss and mass exchange during the common envelope 
phase (CE), when the most massive star, already a neutron star 
or a black hole accretes mass from its companion, a red giant. 
The delay is the sum of the evolution time between the birth 
of the two massive progenitors and the formation of the com- 
pact binary (t ew i), and the merging time, or the time it takes 
for the two stars to coale sce through the emission of GWs (t mr ) 
dPeters & Mafhewsll 1963b . 



(i) 



In order to model the population of compact object bi- 
naries we used the S tarTrack population synthesis code 
dBelczvnski et al. 2002) to perform a suite of Monte Carlo sim- 
ulations of the stellar evolution of stars in environments of two 
typical metallicities: Z = Z Q = 0.02 and Z = 10% Z G = 0.002 
jBelczvnski et al.|[2010bl) . In these calculations we emplo y the 
recent estimates of mass loss rates (Belczvn ski et al.l2010al) . We 
calculate a population of 2 million massive binary stars, track- 
ing the ensuing formation of relativistic binary compact objects: 
double neutron stars (BNS), double black hole binaries (BBH), 
and mixed systems (NSBH). Our modeling utilizes updated stel- 
lar and binary physics, including res ults from supernov a simula- 
tions and compact object formation (iFrver et al.l l2012). incorpo- 
rating elaborate mechan isms for treating stella r interactions like 
mass transfer episodes (IBelczvnski et al.f2008l) or tidal synchro- 
nization and circularization dHutll 19811) . We pu t special empha - 
sis on the common envelope evolution phase (Webbink 1984), 
which is crucial for close double compact object formation as 
the attendant mass transfer allows for efficient hardening of the 
binary. 

An important parameter describing the common envelope is 
the A coefficient. It is a measure of how strongly the donor's en- 
velope is bound to the core. We use the realistic cal culations 
of the common envelope coefficient A performed by IXu & Lil 
(120101) . which now depends on the evolutionary stage of the 
donor, its mass at Zero Age Main Sequence (ZAMS), the mass 
of its envelope, and its radius. The orbital contraction occur- 
ring during the common envelope can be sufficiently efficient 
to cause the individual stars in the binary to coalesce and form 
a single highly rotating object, thereby aborting further binary 
evolution and preventing the formation of a double compact 
object. Because of significant radial expansion, stars crossing 
the Hertzsprung gap (HG) very frequently initiate a common 
envelope phase. HG stars do not have a clear entr opy jump 
at the core-envelope transition (llvanova & T aam 2004); if such 
a star overflows its Roche lobe and initiates a common enve- 
lope phase, the inspiral is expected to lead to a coalescence 
(Taam& Sandquist 2000). In particular, it has been estimated 
that for a solar metallicity environment (e.g., our Galaxy), prop- 
erly accounting for the HG gap may lead to a reduction in the 
merger rates of BBH b inaries by ~ 2 - 3 orders of magnitude 
(IBelczvnski et a l. 2007). In contrast, in low metallicity environ- 
ment t his suppression is much l ess severe (~ 1 order of mag- 
nitude; Belczv nski et alJ (l2010bl) ). The details of the common 
envelope phase are not yet fully understood, and thus in what 
follows we consider two set of models, one which does not 



take into account the suppression (optimistic models: marked 
with A), and set that assumes the maximum suppression (pes- 
simistic models: marked with B). In case of NSs we adopt natal 
kick distributions f rom observations of single Galactic pulsars 
dHobbs et alj|2005l) with cr = 265 km/s. However for BHs we 
draw kicks from the same distribution but with the lower mag- 
nitude: inverse propo rtional to the amount of fall back expected 
at the BH formation dFrver & Kalogerall2001l) . In particular, for 
most massive BHs that form with the full fall back (the direct BH 
formation) the magnitude of natal kick is zero (m ore details on 
modeling can be found in iDominik et alJ d2012i) ). Additionally 
we test one more set of models in which the magnitude of NS 
kicks is decreased by factor of 2: cr = 132.5 km/s, as some 
observations and empirically based arguments seem to indicate 
that natal kicks in close binaries are low er than for single stars 
(iDessart et al.1120061: iKitaura et al1l2006l) . The BH kicks are de- 
creased in the similar fashion as in models with the full NS kicks. 
The standard value of cr parameter is denoted by K and smaller 
value by k. The detailed list of models considered in this paper 
is presented in Table [1] 

Detailed properties such as average total mass or average 
chirp mass are listed in Table [2] 

Table 1. The list of models of stellar evolution used in the pa- 
per. Sigma is a parameter of the natal kick distribution, while 
HG specifies whether we have optimistic (+) or pessimistic (-) 
model. For more details see text. 



Model 


Metallicity 


cr [km s 1 ] 


HG 


AZK 


Zo 


265.0 


+ 


BZK 


/ 


265.0 




AZk 


z. 


132.5 


+ 


BZk 


z 


132.5 




AzK 


10% z Q 


265.0 


+ 


BzK 


10% z 


265.0 




Azk 


10% z Q 


132.5 


+ 


Bzk 


10% Z Q 


132.5 





3. The cosmic coalescence rate 

We assume that for each category of binaries (£=BNS, NSBH, 
BBH), the coalescence rate tracks the star formation rate (SFR) 
albeit with some delay td from formation of the massive binary 
to final merger : 

_ A k f 



Zf 



-P(t d ) dtd, 



(2) 



where p t is the star formation rate in M Mpc yr , A is the 
mass fraction of the progenitors in M a ' (see Appendix for more 
details), P k {td) is the probability distribution of the delay td, Zf 
is the redshift at which the progenitor binary forms, and z is the 
redshift at which the two compact objects coalesce. The factor 
(1 + Z/) _1 converts the rate in the source frame into a rate in the 
observer frame. The redshifts zj and z are related by the delay 
time td which is the difference in lookback times between zj and 
z: 

dz' 

(3) 



td 



Ho X 



(l+z')E(£l,z') 
The coalescence rate per redshift bin is then given by: 

dR k 



dz 



(z) 



dV 

: Pc(z) — (Z) . 

dz 



(4) 



2 



I. Kowalska et al.: Effect of metallicity on the gravitational- wave signal 



Table 2. Statistical properties of compact binaries used in the 
simulations. For each model we listed average total mass of a 
binary, average chirp mass, and average frequency at last stable 
orbit. 



Model 


< M M > [M ] 


< M c Mrp > [Mo] 


< f ho > [Hz] 


BNS 


BZK 


2.43 


1.05 


1811.26 


BZk 


2.43 


1.05 


1809.39 


BzK 


2.49 


1.08 


1768.23 


Bzk 


2.49 


1.08 


1770.95 


AZK 


2.43 


1.05 


1809.71 


AZk 


2.43 


1.06 


1809.53 


AzK 


2.51 


1.09 


1756.74 


Azk 


2.49 


1.08 


1765.63 


NSBH 


BZK 


9.85 


3.13 


448.39 


BZk 


9.49 


3.04 


469.29 


BzK 


12.45 


3.21 


371.55 


Bzk 


11.77 


3.10 


390.27 


AZK 


9.91 


3.17 


444.89 


AZk 


9.85 


3.13 


449.20 


AzK 


11.66 


3.17 


398.82 


Azk 


11.47 


3.11 


404.41 


BBH 


BZK 


15.60 


6.76 


282.59 


BZk 


15.49 


6.72 


285.12 


BzK 


22.41 


9.54 


215.14 


Bzk 


22.26 


9.50 


214.34 


AZK 


15.56 


6.74 


283.66 


AZk 


15.33 


6.64 


288.50 


AzK 


30.31 


13.08 


188.09 


Azk 


28.51 


12.29 


198.89 



The co-moving volume element : 

d JL (z) = 4n JL«EL, 

dz H E(Cl,z) 
where 

M = — f dz ' 
KZ) Ho Jo £(az')' 

and 

£(Az) = \fo A + Q m (l + z) 3 . 



(5) 



(6) 



(7) 



In these simulations, we used the cos mological model de- 
rived from 7 years of WMAP observations (IKomatsu et al .1201 lh 
and assumed a flat Universe with Q,„ = 0.27, Qa = 0.73 
and Hubbl e parameter Hp = 70.3 k m s 'Mpc -1 . Following 



Regimbaul (1201 lh: IZhu et alJ (l201lh . we used the SFR of 
Hop kins & Beacoml d2006l) T"which is derived from measure- 
ments of the galaxy luminosity function in the UV (SDSS, 
GALEX, COMB017) and FIR wavelengths (Spitzer Space 
Telescope), and is normalized by the Super Kamiokande limit 
on the electron anti-neutrino flux from past core-collapse super- 
novas. This model is expected to be quite accurate up to z ~ 2, 
with very tight constraints at redshifts z < 1 (to within 30-50%). 
In order to account for uncertain ties in the SFR, we also con - 
sidered the model presented by Springel & Hernquist] (120031) . 

However, as we'll see in the next section, and as already noticed 

i 1 j 1 

by Regimbau (201 1), the uncertainties on the star formation his- 
tory affecting essentially redshift z > 2 where sources don't con- 
tribute significantly to the integrated signal, our final results are 
not expected to vary by more than a factor 2, depending on the 



choice of the SFR. The final rate is just the sum of the rates from 
the different categories : 



dR dR BNS dR mBH dR BBH 

-riz) = — -, — (z) + — -, — (z) + — — (z). 
dz dz dz dz 



(8) 



Table 3. Total coalescence rate of all compact binaries per year 
is presented in second column. In the last tree columns we show 
contribution of each type of compact binaries to the total coales- 
cence rate. 



Model 


rate [yr l ] 


rate BNS [%] 


rate NSBH [%] 


rate BBH [%] 


BZK 


154 929 


84.78 


2.09 


13.13 


BZk 


157 559 


86.28 


0.94 


12.78 


BzK 


319 304 


10.70 


11.50 


77.80 


Bzk 


364 564 


7.59 


15.90 


76.51 


AZK 


622 572 


71.62 


3.70 


24.69 


AZk 


784 762 


74.30 


2.54 


23.16 


AzK 


1 606 240 


10.62 


3.56 


85.82 


Azk 


1 842 390 


9.73 


4.55 


85.72 



4. The GW background 

The superposition of the GW signal from sources at all redshifts 
create a background, whose spectrum is usuall y characterized by 
the d imensionless energy density parameter dAllen & Romanol 
119991) : 



1 dp gw 
Pc d\nf 



(9) 



where p gw is the gravitational energy density and p c 
is the critical energy density needed to make the Universe flat 
today. The GW spectrum from the population of extra-galactic 
binaries is given by the expression: 



<VCfl = —fF(f). 

PcC 



(10) 



where F(f) is the total flux and / is observed frequency. In this 
paper, we estimate the total flux (in erg Hz -1 ) as the sum of the 
individual contributions, whose redshift and masses are given by 
the StarTrack simulations: 



F(f) = T 



N 



1 



dE 



{ri^d 2 L (z k ) df 



gH '(/,M\M^), 



(ID 



dE., 



where -jj-( f, M , M c ,z ) is the spectral energy density (in erg 

Hz -1 s~') averaged over orientation for a binary at redshift z k , 
with total mass M k and chirp mass M k . N is the number of 
sources produced by StarTrack and corresponds to the num- 
ber of coalescences in T = 1 yr. The normalization factor T~ l 
assures that the flux has the correct dimension. The evolution of 
coalescing binaries is comprised of three phases, the inspiral, the 
merger and the ring-down of the final black hole. The early inspi- 
ral phase and the ring-down can be modeled accurately by post 
Newtonian expansion and perturbation theory, while our knowl- 
edge of the late inspiral and the merger rely on numerical rela- 
tivity. 

In this work, we follow IZhu et al ] (120111 120 12h and use: 1) 
Taylor T4 templates with 3.0 post Newtonian corrections for 
the population of BNS. The merger and ring-down carries quite 
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a lot of uncertainty due to lack of precise knowledge of the 
amount o f ene rgy emitted in these phases. In a recent paper, 
IZhu et all d2012l) have calculated the contribution of an addi- 
tional Gaussian post merger emission centered at 2 kHz and con- 
cluded that it only affected the high frequency part of the back- 
ground whose contribution to the signal-to-noise ratio (SNR) 
is neg ligible. 2) the phenomenological templates of Aiit h et all 
(2008) including inspiral, merger and ring-down, constructed by 
combining analytical and numerical relativity results in the case 
of non spinning BBH 0, for the population of BBH but also for 
NSBH dShibata et al.ld2009l) have shown that NSBH signals with 
a mass ratio larger than 6 look very similar to BBH waveforms, 
because no tidal disruption happens outside the last stable orbit.) 

In the early inspiral phase, the spectral energy density can be 
described by the quadrupole approximation (assuming circular 
orbit) Q 



,x 10 



dE„ 



df 



■(f,flso,M c ,Z) 



(Gtt)V\M c (1 + z)) 5/3 ! 



1/3 



for /<//«,/(! + z), 



where f\ so = c 3 /(6 VottGM) is the GW frequency at the last 
stable orbit (LSO). 

For each category of binary, we can use the following ap- 
proximate analytical expression, which involves the rates de- 
rived in the previous section and average quantities : 



dE ew - i - j dR j 
4ndKz) ^ (/ ' f ,z)—(z)dz, 



(12) 



where ; = BNS, NSBH or BBH, f,J is the average frequency at 

LSO and M c 1 the average chirp mass of the population. Inserting 
Eq.Q~2]in Eq.[10] we obtain: 



with 



'-sup 



(/) 



^i ( ^ )5/ 3 /2/ 3 

rdupD dR> M dz j„ 

Jo dz ^rfeXl+z)" 3 Z ' 



Zmax if/ < flso /(I + Zmax) 

(//.«///) ~ 1 otherwise 



(13) 



(14) 



This expression provides a very good estimate of the background 
in the low frequency region (before the LSO), and in the case the 
chirp mass and the delay distributions show little correlation. 

5. Results and discussion 

We present the dependence of compact object binary coales- 
cence rate on the redshift in Figure [I] The shape of the distri- 
butions is very similar for all the models, as the distribution of 
the delay does not show noticeable differences from one model 
to the other, expect from differences in the minimal value of the 
delay, which is much smaller than the Hubble time and doesn't 
affect significantly the final distribution. The coalescence rate 
peaks at around z ~ 1.7, and is then shifted toward lower redshift 



2 Templates ex ist also in the case of non-precessing spinning BHs 
(Aiith lit al.l201 1[). Using these templates would have very small effects 
on the results IZhu et al.l d201 lh and wouldn't change our conclusions. 

3 In our calculations we used the formula of iPeters & Mathews] 
d 19631) for non zero eccentricity (the initial eccentricity is provided by 
StarTrack ), but it has no effect in the frequency band of terrestrial 
detectors where most of the orbits have been circularized 
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Fig. 1. Rate of compact binary coalescences as a function of 
the redshift for the different values of the metallicity and for the 
different models Table [T] AZK (gray continuous line with trian- 
gle symbol), AZk (gray dashed line with triangle symbol), BZK 
(gray continuous line with square symbol), BZk (gray dashed 
line with square symbol), AzK (black continuous line with cross 
symbol), Azk (black dashed line with cross symbol), BzK (black 
continuous line with circle symbol), Bzk (black dashed line with 
circle symbol). Black lines correspond to low metallicity and 
gray lines to high (solar) metallicity, continuous lines to low kick 
velocity and dashed lines to high kick velocity. 



compared to the star formation rate whose maximal is rather at 
z ~ 2. The normalization, on the other hand, changes, as it is 
linked to the compact object formation rate specific for each 
model. Typically models A (optimistic) correspond to higher 
rates that models B (pessimistic), and models with low metallic- 
ity lead to higher overall rates than the models with high metal- 
licity. The value of the kick velocity has very small influence 
on the coalescence rate. We must note that the models presented 
are calculated for a few particulars population synthesis mod- 
els. Thus the absolute normalization of the results is uncertain, 
just as any rate calculation in the literature. A detailed discussion 
of the expected rates based on the StarTrack population syn- 
thesis models have been published separately bv lDominiket al.l 
J20T1 . 

The resulting gravitational wave background spectrum for 
the models BZk and Bzk are presented in Figure [3]with the indi- 
vidual contributions of each type of binaries shown separately. 

For BBH and NSBH, Q. gK increases as f 2 ^ (inspiral phase), 
then as f 5 ^ 3 (merger phase), before it reaches a maximum and 
decreases dramatically. The slope change and the peak occur at 
frequencies corresponding roughly to the LSO and the end of the 
ring-down phase at z ~ 1.5 where the coalescence rate is maxi- 
mal. For BNS, the slope drops from 2/3 to 1 /3 at around ~ 100 
Hz, as predicted by the TaylorT4 formula. Adding a Gaussian 
post merger emission would enhance th e BNS contribution be- 
tween ~ 0.5 - 2.5 kHz (see IZhu et ail d2012l) ) but would have 
no effect on the detectability of the background as we will see 
in the next section. The amplitude of the background scales up 
with the rate and the average chirp mass (see Eq.lLTt and is larger 
for the BBH background than for the NSBH and BNS contribu- 
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Fig. 2. Energy density parameter of the background produced by the coalescence of compact binaries for the different values of 
the metallicity and for the different models considered in this paper: AZK (gray continuous line with triangle symbol), AZk (gray 
dashed line with triangle symbol), BZK (gray continuous line with square symbol), BZk (gray dashed line with square symbol), 
AzK (black continuous line with cross symbol), Azk (black dashed line with cross symbol), BzK (black continuous line with 
circle symbol), Bzk (black dashed line with circle symbol). Black lines correspond to low metallicity and gray lines to high (solar) 
metallicity, continuous lines to low kick velocity and dashed lines to high kick velocity. The sensitivity of ET-B, ET-D and Advanced 
LIGO/Virgo are also shown for comparison purpose. 



tions by 1 - 3 orders of magnitude. Consequently, at low frequen- 
cies, the background is dominated by the signal from BBH. After 
the BBH peak, the spectrum decreases dramatically, and we re- 
cover the BNS contribution at high frequencies. The population 
of NSBH contributes very little to the total background since it 
is masked by the BBH signal at low frequencies, and, even in 
the low metallicity case BzK, where it is above the BNS con- 
tribution, its maximum is buried in the decrease after the BBH 
peak. 

The spectra for the different models of TableQ]are presented 
together in Figure [2] All the models show the same characteris- 
tic structure described previously, with a peak between 0.1- 0.6 
kHz, corresponding to the population of BBH, followed by the 
high frequency BNS contribution. The first thing to notice is that 



the kicks do not affect the shape of the GW background spectrum 
significantly. This can be expected as the kicks affect neither the 
chirp mass distribution nor the delay time distribution. The effect 
of the choice of the model A or B (with different treatment of the 
CE with HG donor) affects mainly the normalization of the spec- 
trum. The change of the metallicity of the underlying population 
changes both the level and the shape of the spectrum. The lower 
metallicity population contains more high mass BHs and there- 
fore the low frequency peak shifts to the lower frequencies. The 
position of the high frequency peak remains unchanged, as the 
NS masses do not depend on the stellar metallicity. The level of 
the spectrum shifts up with decreasing metallicity. This is due 
to increase of the formation rate of compact object binaries with 
the decreasing metallicity. 
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Model BZk 



Model Bzk 




10 
f(Hz) 

Fig. 3. Energy density parameter of the background produced by the c oalescence of the different type of binaries : two neutron stars 
(blue continuous line with circle symbol), a neutron star and a black hole (red continuous line with square symbol) two black holes 
(black continuous line with triangle symbol), for model B with low kick and with metallicities Z = Z Q (top) and Z = 10% Z (bottom) 
. The black line with no symbole corresponds to the sum of all the contributions. 



The gravitational wave background as presented in Figure |2] 
is maximal in the frequency band of terrestrial interferometers 
(1 Hz- 10 kHz), and could be accessible by the next generations 
of detectors. The standard detection method of cross correlation 
between two detectors gives a signal-to-noise ratio (SNR): 



SNR = V2T 
An 2 



r 

Jo 



df 



f 6 Pi(f)P 2 (f) 



1/2 



where y is the overlap reduction function characterizing the loss 
in sensitivity due to the separation and the relative orientation 
of the detectors, normalized so that it is one for L-shaped co- 
aligned and co-located detectors, P\ and P2 are the strain noise 
power spectral densities of the two detectors. Notice that this 
method is optimal for a Gaussian stochastic background, which 
is not the case of the backgrounds calculated in this paper (the 
sources do not overlap in the frequency range of t errestrial detec- 
tors). However, the ET Mock Data Challenges (Regimba u et al.l 
120121) has proven that such non Gaussian backgrounds can still 
be detected with no loss of sensitivity by the cross correlation 
statistic. 

In this paper we consider 1) a pair of co-aligned and co- 
located Advanced (second generation detectors such as the 
LIGO Hanford pair in the US. A similar sensitivity can be ob- 
tained by combining the future network of detector pairs world- 
wide (LIGO Hanford and LIGO Livingston in the US, Virgo in 
Italy, GEO in Germany, IndlGO in India, KAGRA in Japan), the 
number of pairs balancing the loss in sensit i vity d ue to their sep- 
aration and relative orientation (IZhu et all (120111) . Mandic pri- 
vate communication). 2) a pair of V-shaped detectors (an an- 
gle of n/3 between the two arms) separated by an angle of 
2n/3, with the sensitivity of the third generatio n triangle detec- 
tor ET. We use two models for the noise, ET-B (iHild et al.f2 008) 
which is about ten times better than the sensitivity of Advanced 
LIGO/Virg o with a low frequency cutoff at 1 - 3 Hz, and ET-D 
dHild et alj|201lb with a significant improvement at frequencies 
below 30 Hz. 



(15) H- 




f (Hz) 



Fig. 4. Contribution to the SNR of frequencies < /, for co- 
aligned and co-located Advanced detectors, separated Advanced 
detectors (the LIGO Hanford/Livingstonpair), and the third gen- 
eration Einstein Telescope (ET-B and ET-D) 



In Fig.|H we show the cumulative SNR as a function of the 
frequency for Advanced and ET detectors for the model BZk 
(considering other models or a simple power law Q gw ~ f 2 ^ do 
not change the results significantly). As one can see, most of the 
SNR comes from the power-law low frequency region, where 
the BBH contribution dominates. 

The SNR is reported in Table |4]fbr all the models, assum- 
ing an integration time of 1 yr. Models with low metallicities 
should be detected by Advanced detectors with SNR ~ 7 - 57, 
so about 2-20 times larger than the detection threshold of 
SNR = 2.5 (for a false alar m rate and false dismissal rate of 
10% dAllen & Romanoll9 99)). Models with high metallicity are 
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accessible at the detection threshold for model A but are out of 
reach for model B by a factor of a few. With Einstein Telescope, 
all the models should be detectable with very high signal-to- 
noise ratio, between SNR~ 39-4155 forET-B and 102-14485 
for ET-D. 



Table 4. For the different models presented in Table [T| signal- 
to-noise ratio obtained by cross-correlating two co-aligned and 
co-located Advanced detectors or a pair of two ET detectors with 
the sensitivity ET-B or ET-D. 



Model 


Ad CC 


ET-B 


ET-D 


BZK 


0.6 


41 


107 


BZk 


0.5 


39 


102 


BzK 


7.9 


568 


1567 


Bzk 


6.8 


490 


1350 


AZK 


3.2 


234 


614 


AZk 


2.7 


195 


512 


AzK 


57 


4155 


14485 


Azk 


52 


3750 


13466 



Another promising detection method that have been tested 
in the context of the ET Mock Data Challenge (Regimba u et alj 
1201 2l) for the population of BNS is to look at the residual of the 
null stream which contains no GW signal and arises from the 
sum of the three ET detectors. 



6. Conclusions 

In this work, we investigated the effect of the metallicity on the 
background from extra-galactic compact binary coalescences. 
We used the Star/Track population synthesis code to perform 
a suite of Monte Carlo simulations of the stellar evolution of 
binary systems in environments of two typical metallicities: 
Z = Z Q = 0.02 and Z = 10% Z = 0.002. We found that for 
all models, the background is dominated by the binary black- 
hole contribution in the frequency range where terrestrial detec- 
tors are the most sensitive, with a common structure : an evolu- 
tion as £2 gw ~ f 2 ^ corresponding to the inspiral phase, then as 
Sl gw ~ f 5 ^ corresponding to the merger phase, a maximum at 
0.6- 1 kHz, and finally at high frequencies, the final stages of the 
binary neutron star contribution in the post Newtonian regime 
with Cl gw ~ / I/3 . A smaU metallicity Z = 10% Z Q = 0.002 pro- 
duces more massive systems leading to a larger amplitude and a 
shift of the peak frequency to lower values. 

All our models are detectable by Einstein Telescope after a 
year of observation with a very high signal-to-noise ratio (be- 
tween 40 - 14000). Models with low metallicity are also de- 
tectable by Advanced LIGO/Virgo, while models with solar 
metallicity are close to the detection threshold with an SNR 
betwe en 0.5 - 3., in a g reement w i th the estimate of detectabil - 
itvbvlWuetal.1 (120121) : IZhu et alJ (120111) : iMarassi et al.1 (1201 lh : 
Zhu et al. (2012). Since the background from compact binary co- 
alescences has good chance to be detected in the near future with 
the second generation of detectors, having realistic estimates of 
its shape and amplitude is crucial in order to optimize the search 
and later interpret the results. The LIGO/Virgo stochastic group 
and the ET Science Team have started Mock Data Challenges 
based on the models presented in this paper whose results will 
be presented in future papers. 
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Appendix A: Normalization of the cosmic 
coalescence rate 

To normalize coalescence rate, we need to take into account var- 
ious factors. First of all, our simulations can not be treated as a 
complete sample of realistic stars. We do not cover whole pa- 
rameter space, in particular the mass range is narrower than in 
reality. Also, not all stars that are born in binaries will end up as 
a compact binaries. Most of them will not survive the common 
envelope stage. 

Coalescence rate for one particular binary i is given by fol- 
lowing formula: 

Pcj(z) = —j—n bin f s , (A.l) 
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where N S i m is the total number of binaries that started their 
evolution on ZAMS in the Star/Track code. In our case that 
number is 2 x 10 6 , n b j„ is the number density of stars in binaries, 
f s = is the fraction of simulated binaries. 

In the Universe not only binares are formed. Total Star 
Formation Rate (p*) is a sum of a contribution from binaries and 
single stars. 



p»=p« bin +p, sin [M e Mpc yr ]. 
Number density will be defined as: 



n k = 



P*k 



<M k > 



[Mpc yr ], 



(A.2) 



(A.3) 



where k="bin" or "sin" and < M k > is the average mass of a stat- 
in a binary, if k — bin, or a single star, if k — sin. 

Another important parameter is binary fraction, which is de- 
fined as a fraction between number of binary systems and total 
number of all systems (binary and single). 



ft 



Nbin 



0.5n bil 



Nbin + N si „ 0.5n hin + n s 



(A.4) 



where N b i n is the number of binaries, not the stars, so we 
need factor of 0.5 in the second par of the equation. We set binary 
fraction to be f b = 0.5, as it is for a current stellar populations. 

By transforming Eq. IA.4l we obtain number density of a sin- 
gle stars: 



- 1 \zA 

tlsin — Whin ~ ) 
1 Jb 



(A.5) 



then, combing eq. IA.2I and IA.3I we have number density of 
stars in binaries: 



and M m 



M max Then, the mass of the secondary is deter- 



mined by mass ratio, which is taken from a flat distribution. Mass 
of the second component has to be greater than 3M . Fraction of 
simulated stars can be express by: 



£~ dAW(M) £ dqO(q) 



(A. 12) 



where q min (M) = -|. 

Normalization A* (k=BNS, NSBH or BBH) can be define as: 



1 V vim 



(A.13) 



Af is the number of k-type binaries in our simulations. 
Detailed list of that values for all models is presented in Table 
IAT1 



Table A.l. Number of binaries of each type produced in all mod- 
els. 



Model 


N BNS 




N tiHU 


BZK 


1618 


43 


395 


BZk 


1994 


56 


441 


BzK 


509 


475 


2873 


Bzk 


464 


885 


3175 


AZK 


4940 


326 


1721 


AZk 


6710 


404 


2164 


AzK 


1698 


736 


15312 


Azk 


1769 


1283 


17613 



ribi„ = p. t-7 . ( A.6) 

< M sin > +2 < M bin > 

(A.7) 



bm < M sin > +2 < M bi „ > 

As a next step we need to compute average masses of single 
star and of those that are components of a bi nary. We used I nitial 
Mass Function (IMF) proposed by Salpeter (Sal peterlll955l) and 
flat mass ratio distribution. 



^(M) = KM- 2 \ (A.8) 
q>(q) = 2q, (A.9) 

where K~ l = J^""" M- 2J dM, M mi „ = O.O8M , M max = 
150.1/ . 

Using IMF we can compute average masses of a single stars 
and components of a binaries. 

pM,„„ 

<M si „>= I M x ¥(M)dM, (A. 10) 

Jm„,„ 

<M bi n>=- I M x V(M)dM I dq(\ + q)<5>(q), (A. 11) 
2 Jm w „, Jo 

The form of f s is a consequence of a specific way of drawing 
masses of the components in StarTrack. First mass of the pri- 
mary is drawing from Salpeter distribution with M m ,„ s , m = 5M© 
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